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We study the spectrum and eigenmodes of the QCD Dirac operator in a gauge background given 
by an Instanton Liquid Model (ILM) at temperatures around the chiral phase transition. For two 
massless quark flavors we observe that at the chiral phase transition, both the low lying eigen- 
modes and the spectrum of the QCD Dirac operator resemble those of a disordered conductor at 
the Anderson transition. In particular we present results showing that the eigenvectors are mul- 
tifractal and the spectral correlations are well described by critical statistics. We also find an 
additional temperature dependent mobility edge located in the bulk of the spectrum with proper- 
ties quantitatively similar to those of a 3D disordered conductor at the Anderson transition. This 
feature is present in both quenched and unquenched simulations. Finally we argue that our find- 
ings are not in principle restricted to the ILM approximation and may also be found in lattice 
simulations. 
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1. Introduction 

The spontaneous chiral symmetry breaking (S;t^SB) in QCD and its eventual restoration at 
finite temperature is one of the most intensively investigated topics in hadronic physics. QCD 
models whose gauge configurations are given by an interacting liquid of instantons ^ provide an 
adequate theoretical framework to understand S;^SB. Thus, based on the semiclassical picture of 
a QCD vacuum dominated by instantons, it has been suggested |Q] that the S;fSB in QCD and 
the phenomenon of conductivity may have similar physical origins. Conductivity in a disordered 
sample is produced by electrons that although initially bound to an impurity may become delocal- 
ized by orbital overlapping with nearby impurities. Similarly, in the QCD vacuum, the zero modes 
of the Dirac operator though initially bound to an instanton may get delocalized due to the strong 
overlap with other instantons. As a consequence, the chiral condensate becomes nonzero and chiral 
symmetry is spontaneously broken. In this paper we provide evidence that these analogies can be 
extended to describe the chiral phase transition in QCD at finite temperature. Specifically, in the 
context of an Instanton Liquid Model (ILM), we show that around the temperature that the chiral 
phase transition occurs the low lying eigenmodes of the QCD Dirac operator undergo an Ander- 
son transition (AT), namely, a localization-delocalization transition, characterized by multifractal 
eigenstates and critical statistics ||3|j. 

For the sake of clarity we briefly summarize the properties of a disordered system at an AT. 
In three and higher dimensions there exists a mobility edge separating localized from delocalized 
states. In the delocalized region, eigenfunctions are extended through the sample and the level 
statistics are described by random matrix theory (RMT). In the opposite limit, eigenfunctions are 
exponentially localized and the spectral correlations are described by Poisson statistics. Around 
the mobility edge, the region where the AT occurs, the eigenstates are multifractal, meaning that 
moments of the eigenfunctions exhibit anomalous scaling with respect to the sample size iQ]. Typ- 
ical features of the level statistics at the AT, dubbed critical statistics [|5|], include scale invariant 
spectrum [|6|], level repulsion and asymptotically linear number variance. Thus the level spacing 
distribution goes to zero as the spacing goes to zero as in a metal but long range correlators such 
as the number variance T}{€) = {{Ni — {Ne))^) ~ for ^ ^ 1 {N( is the number of eigenvalues 
in an interval of length £) are asymptotically linear, as for an insulator (x = 1), but with a slope 
X <^ 0.27 for a 3D Anderson model [Q]). Below we show under what circumstances the above 
features are also found in the ILM at finite temperature. 

2. The instanton liquid partition function: Technical details 

Our starting point is the Euclidean QCD partition function, 

Zinst = det^' (0 + m)e-^™ , (2.1) 

where the integral is over the collective coordinates of the instantons, the Dirac operator is denoted 
by p and the Yang-Mills action is denoted by 5ym- For each instanton we have 12 collective 
coordinates (for three colors). The fermion determinant is evaluated in the space of the fermionic 
zero modes of the instantons (for further details we refer to the review [^). For simplicity we have 
kept the density fixed atN/V = I fm^^ as the temperature is increased. This is justified since it was 
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Figure 1: Level spacing distribution P{s) in the 
bulk for the quenched ILM at T = 200 MeV for 
different volumes. The inset shows \ogP{s) for 
the tail of the data. The best fit (solid hne) corre- 
sponds to a slope — 1 .64. 




Figure 2: Spectral rigidity in the bulk for the 
quenched ILM at T = 200 MeV for different vol- 
umes. The result has very little size dependence 
and agrees well with the prediction of critical 
statistics (solid line). 



already observed that the density of instantons at the transition is still sizable, around 0.6 fm^^, and 
the transition is not expected to be driven by the decrease in density. We are therefore using this as 
a qualitative model for the QCD chiral phase transition even though the transition temperature will 
differ from the lattice prediction. 

We report results for Nf = (quenched) and Nf = 2, m„ = = (unquenched). The effect 
of a nonzero mass or a strange quark will be discussed elsewhere [^]. The partition function (2.1) 
is evaluated by means of a Metropolis algorithm. In all simulations we performed between 500 
and 1000 equilibration updates followed by 2000 to 10000 measurement updates. In the quenched 
case results are presented for ensembles of up to A/^ = 5000 instantons and anti-ins tantons. In the 
unquenched case, the simulations are slower and we have only investigated ensemble of up to 
A'^ = 504. All the spectral correlators are calculated from the unfolded spectrum. This procedure 
scales the eigenvalues so that the average density on the scale of several level spacings is one. 



3. Localization in the Instanton Liquid Model: The quenched case 

Here we examine the properties of eigenstates in the center of the spectrum (bulk). Since only 
the smallest eigenmodes are responsible for S;^SB, the mobility edge seen here is not related to the 
chiral transition (strictly speaking there is not a chiral transition in the quenched case anyway). We 
instead study it because it is a very clear example of a mobility edge in the ILM and we find it has 
similar features to the mobility edge near origin, but since it affects a larger range of eigenmodes 
we can get better statistics. By using the finite scaling method introduced in |^ we have observed 
a mobility edge separating localized from delocalized eigenstates in the range T ~ 150 - 250 MeV. 
As the temperature decreases the location of the mobility edge moves to the end of the spectrum. 
For r < 170 MeV the results are less reliable since the mobility edge is located almost in the end 
of the spectrum where truncation effects are larger. Following the literature in disordered systems 
we have investigated the temperature T 200 MeV such that the mobility edge is located around 
the center of the spectrum. The level statistics in this spectral region have the signatures of an 
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Figure 3: Level spacing distribution P{s) for the quenched ILM at = 5000 and T = 200 MeV for different 
spectral windows. The ranges indicate the eigenvalue numbers plotted with being the smallest eigenvalue. 



AT including a spectrum that is scale invariant to high degree. As observed in figures [T] and ^ 
both short range statistics such as the level spacing distribution, P{s), and long range correlators 
such as the spectral rigidity ^^{l) = /o (^^ —2£^x+x^)L^{x)dx do not depend on the system size 
for volumes (number of instantons) ranging from A'^ = 500 to N = 5000. Level repulsion is still 
present as for a metal but the spectral rigidity (see figure ^) is asymptotically linear with a slope 
that corresponds to ;^ ~ 0.29 ±0.02 in fair agreement with the value for a 3D disordered system at 
an AT of X ^ 0-27 ±0.02 Furthermore the exact form of the spectral rigidity follows closely 
the prediction of critical statistics, equation (31) in [10] with h, a free parameter, set to a value of 
0.62. As observed in figure ^ the level statistics outside the critical region are RMT like in the 
region between the origin and the mobility edge and Poisson like between the mobility edge and 
the end of the spectrum. 

We now present the eigenvector properties. A signature of an AT is the multifractality of 
the eigenstates as measured by the moments Pg = J d'^r\Y(j')\^'^ °^ i^-D^[q-\) ^ where Dq is a set of 
different exponents describing the transition and Loc N^/^ is the spatial system size. In order not to 
consider eigenstates outside the critical region we have taken only 10% of the eigenvectors around 
the center of the spectrum for T = 200 MeV. By fitting P2 (also known as the inverse participation 
ratio or IPR) for different volumes, we have obtained a value of D2 ~ 1.5 ±0.1 in agreement with 
that of a 3D disordered system at the AT [11]. 

Our numerical analysis shows that at T = 200 MeV there is mobility edge around the cen- 
tral part of the spectrum in the quenched ILM. The level statistics have all the signatures of an 
Anderson transition: scale invariance, level repulsion and sub-Poisson spectral rigidity; and the 
eigenstates appear to be multifractal. We have also observed similar- critical features in the region 
close to the origin at temperatures around 110-140 MeV. However the analysis is complicated by 
the accumulation of very small eigenvalues that are present in the quenched ILM. More details of 
this region will be given elsewhere [^. Instead we present results near the origin for two massless 
flavors where the accumulation of very small eigemnodes is suppressed and the analysis is cleaner. 
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Figure 4: Chiral condensate and inverse partic- 
ipation ratio (IPR) of the lowest eigenmode for 
the ILM with two massless quarks. 



Figure 5: Level spacing distribution P{s) at the 
origin for the unquenched ILM at T = 115 MeV 
for different system sizes. 



4. Chiral restoration and Anderson transition: The unquenched case 

We now study eigenvalues and eigenvectors of the QCD Dirac operator in an ensemble with 
two massless quark flavors. In this case the chiral condensate associated with S;^SB is finite. 
Around the critical temperature, the condensate approaches zero signaling the chiral phase transi- 
tion. Since we have two massless flavors we expect to find a second order chiral phase transition. 

We focus our attention on the spectral region close to origin since our aim is to relate the chiral 
phase transition with a transition to localization (AT) of the QCD Dirac operator eigenmodes. By 
looking at the chiral condensate versus temperature for a range of system sizes (figure |^ we see 
that the condensate approaches zero around temperatures of 1 15 - 120 MeV in a manner consistent 
with a second order phase transition. At the same time the condensate is falling we see the IPR 
of the lowest eigenmode begin to rise signaling a transition to localization. This finding strongly 
suggests that both phenomenon are intimately related. 

Next we examine the nature of the localization transition. First we locate, by using the finite 
size scaling method [^, the AT for the lowest lying eigenmodes at T ~ 115 - 120 MeV. At this 
temperature we see that eigenvalue statistics such as the level spacing distribution (figure ^) are 
nearly scale invariant, there is level repulsion but the tail of P{s) is exponential as at the AT. We also 
looked at the scaling of Pq in the lowest eigenmode at T = 120 MeV for system sizes ranging from 
=63-252. This yielded a set of fractal dimensions D2 = 1.3 ±0.2, D3 = 0.9 ±0.2, D4 = 0.7 ±0.3 
and D5 = 0.7 ± 0.3 showing that the eigenstates are multifractal. Although higher volumes would 
be desirable to fully confirm the multifractal scaling of the eigenstates, it is encouraging that the 
IPR of comparatively similar temperature {T ~ 100 Mev) is comparable to the one for a fully 
metallic sample. We also mention that we observed a mobility edge on the bulk of the spectrum at 
a slightly higher temperature with properties similar to the quenched case. 

At this point a natural question to address is whether the results here reported concerning 
the Anderson transition in the ILM may also be present in more realistic models of the strong 
interactions. Although no conclusive answer can be given to this question we would like to point 
out the reasons why we have found an Anderson transition in the ILM. As is known, unlike the 
zero temperature case where the decay is power-law, the fermionic zero modes in the field of 
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an instanton at nonzero temperature (usually called a caloron) have an exponential tail e^''^ in 
the spacial directions and are oscillatory in the time direction [||]. This suggests that the overlap 
among different zero modes is essentially restricted to neai^est neighbors in the spacial directions. 
However, in the time direction different zero modes strongly overlap due to the oscillatory character 
of the eigenmodes. This situation strongly resembles a 4D disordered conductor in the tight binding 
approximation (only nearest neighbor hopping) with one dimension (time) much smaller than the 
rest so the system can be considered effectively three dimensional. It is well established that such 
a system may undergo an AT depending on the disorder strength. In our case the disorder role is 
played by temperature since the wavefunction decays as ^ e^''^ . From this discussion it is natural 
to find an Anderson transition in the ILM for a particular value of the temperature. The principal 
ingredient to reach the AT is an exponential decay of the eigenmodes explicitly depending on the 
temperature together with the possibiUty to tune the effective range of the exponential in order 
to reach the transition region. Any theory with these features very likely will undergo an AT for 
some value of the parameters. It is therefore quite possible that the QCD vacuum has similar 
properties. Even if the objects responsible for localization are not the classical instantons (see the 
topology section of these proceedings for other options) this scenario could still be realized if the 
interactions behave in a similar way. We refer to for a more detailed account of this interesting 
issue. 

In summary we have found that in the unquenched ILM the chiral restoration transition occurs 
around the same temperature as the localization transition close to the origin. From a physical 
point of view this is not surprising since the appearance of the condensate is deeply linked with 
the delocalization of the zero modes in the instanton vacuum due to long range instanton induced 
interactions. It will be very interesting to repeat these same studies in lattice simulations where 
the chiral transition also coincides with a deconfinement transition and additionally to see if the 
detailed structure of the topological objects affects the nature of the localization transition. 
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